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Combining infrared and low-temperature asymptotes in Yang-Mills theories 
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We demonstrate that power-like non-perturbative behavior of gluon and ghost propagators in 
infrared limit of Yang-Mills theories can provide at finite temperatures T a negative T 4 -contribution 
to the pressure and energy density. Existence of mass gap implies then new relations between infrared 
critical exponents of gluon and ghost propagators. 
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Zeros of Green's functions play important role in var- 
ious condensed matter systems ranging from superfluid 
3 He to unconventional superconductors Under cer- 
tain circumstances poles of Green's functions in momen- 
tum space - which determine the excitation spectrum of 
the system - can undergo transition to their opposites, 
that is, zeros. A known example of such a system is the 
Mott insulator. In this pole of a one-particle 

Green's function becomes a zero as the system passes 
through the Mott transition Q • In field theory, a similar 
example is provided by strong interactions. As was first 
argued by Gribov in the Landau gauge a perturba- 
tive pole of the gluon propagator is converted into a zero 
at the vanishing momentum 0, [B| . This phenomenon is 
confirmed by both analytical calculations and numerical 
simulations 0. The effect is related to the violation of 
reflection positivity and is a signature of the color con- 
finement. The aim of this Letter is to demonstrate that 
anomalous behavior of the gluon propagator in the in- 
frared region is directly related to the anomalous contri- 
bution of massless degrees of freedom at vanishing tem- 
peratures. The very existence of a mass gap in confining 
theories constitutes then a new kind of analyticity which 
constrains infrared asymptotes of the propagators. 

We consider pressure P and energy density e = E/V 
of strongly interacting Yang-Mills fields at temperature 
T and in volume V. In the limit of vanishing coupling, 
or in the tree approximation: 



e freo (T) = 3P free (T) = N d . f . C SB T 4 , 



(1) 



where Csb=tt 2 /30 is the Stefan-Boltzmann factor. There 
are N±f. — 2(A r ( ? — 1) degrees of freedom corresponding 
to N% — 1 gluons with two transverse polarizations. 

Non-perturbative strong interactions are manifested 
in suppression of the energy density and pressure at 
T <C T c , where T c is the temperature of the confinement - 
deconfinement phase transition: 



In particular, such a suppression is established via lattice 
simulations both in the pure gluonic SU (3) case Q and 
in a more realistic case which incorporates light quarks 
as well [8]. Hereafter we will concentrate on the pure 
gluonic SU(3) gauge system. 

The suppression ([2]) is a consequence of two intimately 
related properties of QCD, namely, confinement of color 
and mass-gap generation. Massless gluons are confined 
into glueballs which are massive. As a result, all the 
thermodynamic quantities are to be suppressed at low 
temperatures as 0(exp{— M/T}), where M is the mass 
of the lightest glueball, for a related discussion see [§] . 

The simplest quantity for a homogeneous system in the 
thermodynamic equilibrium is pressure, 



P= - log £ = P free + Pint 



(3) 



e(T) 

rp4 



< 1 



P(T) 



< 1 for T < T c 



(2) 



where Z is the partition function and we explicitly sep- 
arate the tree level contribution Pf ree given by Eq. |T]). 
The second term P mt is a correction due to the interac- 
tions. The interaction part of the thermodynamic quan- 
tities can be written as a loop expansion, in terms of the 
full propagators and vertices, see, e.g., Ref. [Tpj ] . 

In order to make explicit calculations one generically 
needs a gauge fixing, and below we usually refer to the 
Landau gauge which - being both Lorentz- and color- 
symmetric - is one of the best studied gauges nowa- 
days. However, our approach is not limited to the Landau 
gauge, and, as we will see below, is in fact very general. 
As a consequence, the final result - a relation between the 
infrared exponents of the dressed gluon propagator and 
the dressed ghost propagator should be gauge invariant, 
provided that the infrared exponents (which may vary 
from gauge to gauge) are defined as we discuss below. 

Before going into details we would like to make a gen- 
eral comment about the gauge invariance. Despite of the 
fact that Green's functions in non-Abelian theories are es- 
sentially gauge- variant, their properties may be directly 
linked to the confinement. The best known example of 
this kind is provided by Ref. Q which relates the linear 
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confining potential increase of the color charge interac- 
tion at large distances to cancelation of the infrared pole 
of the gluon Green's function. 

At finite temperatures, the gluon propagator D ab v and 
the ghost propagator G ab in the Landau gauge are pa- 
rameterized by three functions Dj<, D^, and Dq: 



gat 



PL D T {v,Pi) + P± D L (p,p 4 ) , (4) 



G ab (p, P4 ) = -6 ab D G (p, Pi ). 



(5) 



where P^ v and Ph, are projectors onto spatially trans- 
verse and spatially longitudinal parts [ill ]. 

A conspicuous feature of the propagators is the anoma- 
lous dressing of the gluon and ghost propagators in the 
infrared region [6]. In the low temperature limit: 



D t tf) = (p 2 f> ^(p*) as p 2 







(6) 



where 7$ are infrared exponents and the functions Hi are 
finite at p 2 — ► in all the cases, i = T,L, G. 

The infrared behavior of the propagators is best stud- 
ied at zero temperature. First of all, the gluon infrared 
exponents are degenerate, jd — 7t = 7l- We have 
already mentioned the Gribov's scenario for 70 > 
[H, EL 01- Moreover, according to the Kugo-Ojima crite- 
rion of confinement [12j the ghost propagator is enhanced 
in the infrared, compared to its perturbative value: 



7£> + 27c = . 



(7) 



There is some support for the validity of this rela- 
tion, both from the analytical studies of the Dyson- 
Schwinger equations @, and from numerical, lattice 
data [14j, ll5( . However, final conclusion on the validity 
of Eq. ((7]) seems to be not yet reached, see in partic- 
ular Refs. HI, 3, 17, A cautionary remark comes 
also from a confining Abelian gauge model where one can 
prove analytically [19j that Eq. ([7]) does not hold. 

Turning to non-zero temperatures, let us note that 
the infrared properties of the gluon propagators have al- 
ready been used to study thermodynamics of Yang-Mills 
plasma at T > T c in the Coulomb gauge (20|. The idea 
is that constraining the configuration space through the 
gauge fixing condition leads to a change in the disper- 
sion relations for gluons, or position of the poles of the 
propagator. In the Landau gauge, a relation between 
the non-perturbative Green's functions and the thermo- 
dynamic potential was discussed in Refs. 21, . Below 



we concentrate on zeros rather than poles of the prop- 
agators and demonstrate that zeros give rise to new T 
terms in the equation of state. 

To illustrate the basic idea consider first a toy model 
describing one degree of freedom, or a real scalar field 
4>(x) with a quadratic action: 

S[4>] = -~ f /T dT J d 3 x cf>(x) (x - y) cj>{y) , (8) 



where the propagator V(x — y) is defined as 

V{x -y)= {cf>{x)ct>{y)) , (9) 

and we use the imaginary time formalism. The partition 
function can be calculated in the standard way [111 ]: 



Z = / D^e" 



exp{iTrlog(T 2 P)}, (10) 



where we omit the irrelevant prefactor. Taking the trace 
over all the states, one gets for the pressure ([3]): 



p= T [d^p 



2 J {2irf 



^log[T 2 2?(p,p 4 ) 



, (11) 



where the sum runs over the Matsubara frequencies, 
uj n = 2nnT. The energy density is given by 



e = -T 



d 3 p y, 9l0g[p|2?(p,p4)] 



(2tt 



<9l0gP4 



(12) 



The temperature depende nt p art of this sum can be 
evaluated following Refs. 111. |23||: 



e(T) = 



1 f d 3 p + 7 + , £ , , > F(Po) + F (-go) 
Vi JWf ' dpofTM 5 (13) 



— ico-\-e 



FeiPo) 



d\0g[plV(p,p 4 )} 



d log pi 



(14) 



Pi = -Wo 



where Jt(po) = l/(e p °/ T — 1) is the Bose-Einstein distri- 
bution for a particle with energy to = po- Equation (|13p 
is valid provided the analytical function F(po) does not 
have poles at purely imaginary po- Since the integrand 
in Eq. ()13|) is fast converging as po — * +00, the contour 
of integration can be closed as a semicircle in the right 
half of the po-complex plane, reducing the integral (fl3|) 
to contribution of poles of the function F e (p ). 

Poles of the propagator D become poles of the func- 
tion F s (po) and contribute to e(T). This agrees with our 
intuition since the energy spectrum of the model JHJ) is de- 
termined by the poles of the propagator V. However, the 
crucial point is that not only the poles of the propagator, 
but also its zeros contribute to the energy density (|13[) . 
(fT4|) . In fact, zeros of the propagator V become poles of 
the function F e (j>o) as well. 

Imagine that the propagator f9| in the infrared region 
has the same criticality as the gluon or ghost propaga- 
tors © in Yang-Mills theory: 



D {l) =const(p 2 )7(p 2 + m 2 ). 



Then we get from Eq. (fT4j) the expression 

,2 



F^(po) 



E__ 7 

-Po P 2 Po 



(15) 



(16) 
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with two poles at the Repo > half of the complex p - 
plane: at po — u> p and at po = |p|. The residues, respec- 
tively, are — w p /2 and 7|p|/2. The energy density (H^J) is 



e(T) 



s bcc {T,m) - 7 e frcc (T,0), (17) 



£ free (T,m) = J^^fr^). 



(18) 



Here e flcc is the energy density of free relativistic particles 
of mass to. 

It is clear from Eq. (fT?)) that the infrared suppression 
(with 7 > 0) in the particle propagator acts oppositely 
to the pole: the suppression "subtracts" massless degrees 
of freedom (with the dispersion relation corresponding 
to a relativistic massless "particle" lo p = |p| and with 
£f ree (T, 0) oc T 4 ) while a pole "supplies" them to the spec- 
trum. The same effect is also from the expression for the 
pressure (fTTj) because \ogD^ = j\og(p 2 )—\og(p 2 + m 2 ), 
up to a momentum-independent term. Note that a zero 
of the propagator corresponds to a negative pressure and 
energy density. Thus, vanishing of the propagator in the 
infrared with critical exponent 7 acts as a "reservoir" for 
7 massless degrees of freedom where the quantity 7 is, 
in general, non-integer. As we see from our illustrative 
example, the low-temperature asymptotics of the ther- 
modynamic functions are determined by the infrared be- 
haviour of the particle propagator. 

Let us now turn back to Yang-Mills theory. In the lead- 
ing, one-loop approximation (in terms of dressed propa- 
gators) pressure in Yang-Mills theory at vanishingly low 
temperatures is given by the expression (jTTJ) with 



V{p 2 



D L (p 2 )D 2 T {p 2 )[D L (p 2 ) + Dt(p 2 )} 

d 2 g {p 2 ) 



(19) 



The ghost structure function Dq in the denominator of 
Eq. (fT9)l corresponds to the Faddeev-Popov determinant 
which subtracts unphysical degrees of freedom. Indeed, 
if one neglects the interactions, g = 0, then 



D { [ ee = D£ ee = D% cc = l/p 



-»frcc 



free 



(20) 



and therefore the energy and pressure are given by Eq. (TT|) 
with Nd.f. = 2(N 2 - 1), as expected. 

However, if one takes into account the infrared dressing 
of the propagators ([5]), then one gets for the pressure 



P = (2 - 7i - 27 T + 2 7G - min[7 L , 7T ]) 
■(N 2 - 1)C SB T 4 + . . . , 



(21) 



where ellipsis denote subleading 0(T K ) terms with k > 4, 
associated with cuts in the complex p -plane. The sup- 
pression ^ of thermodynamic quantities at low temper- 
atures implies a relation between the infrared exponents: 



Since the functions Hi, i = L,T, G, are regular at p 2 = 
the poles and/or zeros of these functions do not con- 
tribute to a T 4 term at very low temperatures. Indeed, 
poles and/or zeros at momenta Repo > would corre- 
spond to massive-like contributions which are exponen- 
tially suppressed by the Bose-Einstein factor as T — > 0. 

Coming back to the question of gauge dependence of 
our results, we mention that the form of the propaga- 
tors (J6|) assumes local Lorentz isotropy of the gauge con- 
dition at zero temperature. However, in a general case 
(like, for example, in the Coulomb gauge) or/and at a fi- 
nite temperature (as we discuss in detail below) the form 
of the propagator is not rotationally invariant, and the 
role of the infrared critical exponent is played by the de- 
gree 7 of the zero, corresponding to the free dispersion 
relation pg cc = |p|. In the Coulomb gauge effects of the 
Fundamental Modular region lead to a substantial mod- 
ification of the gluon dispersion relation [2(j. Namely, 
relation between the energy and momentum of gluon 
is no longer linear at low momenta. Thus, the gluon's 
contribution to the equation of state in this gauge does 
not contain T 4 terms in the limit of low temperature, in 
agreement with lattice simulations (20j . 

The relation (|22|) provides a new constraint on the crit- 
ical exponents. Existing lattice data does not allow to 
check it directly. One can try to unify (f2"Tj) with other 
relations suggested in the literature. Assuming validity 
of ([7]) and 7l = 7 t we get a system of constraints which 
allows to uniquely fix all the exponents and get what we 
would call simplified solution, 773 = 2/5 and 7g = —1/5. 

The positivity of 70, exhibited by this solution, is con- 
sistent with the infrared suppression of the gluon propa- 
gator expected on other grounds Q. Also, the ghost ex- 
ponent is consistent with the numerical result of Ref . 17 1 . 



7 l + 277 - 2jg + niin[7L, 7t] -2 = 0. 



(22) 



However, for the simplified solution the gluon infrared ex- 
ponent has too small value to guarantee vanishing of the 
gluon propagator at p = 0. Also, the simplified solution 
is in variance with results of Refs. [HI, EH ■ One of the rea- 
sons could be violation at the presently available lattices 
of the equality = It which is granted at zero tem- 
perature. Indeed, according to numerical results of Ref. 
22] ] the transverse-longitudinal degeneracy is not reached 
yet even at the lowest non-zero temperatures available 
now. One might guess that the discrepancy between the 
simplified solution and the existing simulations is rooted 
in the fact that the lowest currently available non-zero 
temperatures are still too high to be compared with our 
analytical predictions. 

One can ask a question whether the relation ([22|l be- 
tween the critical exponents is affected by the presence 
of quarks. As far as the quarks are massive both numer- 
ical 25 1 and analytical results [24| indicate absence of 
the infrared particularities (zeros or poles) in the quark 
structure functions in the Landau gauge. Consequently 
our relation (|2"2"|) is insensitive to the presence of the mas- 
sive quarks. This conclusion agrees nicely with the re- 
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suit of Ref. 24|. However, for exactly massless quarks 
the T 4 -terms do appear in the thermodynamical quan- 
tities in the T — > limit because of the presence of ex- 
actly massless Nambu-Goldstone particles in the spec- 
trum. Then the right hand side Eq. (f2"2")l should be equal 
to — 7Vng/(-^c — l)i where -/Vng is the number of the 
Nambu-Goldstone bosons. 

To summarize, the zeros in the momentum-space 
Green's function are as important for the thermodynam- 
ics as the poles. Although zeros do not correspond to any 
states in the spectrum they contribute to the equation of 
state. The contribution of zeros in the gluon propagator 
to coefficients in front of T 4 terms in the energy den- 
sity and pressure is negative. This observation can be 
interpreted in such a way that zeros correspond to con- 
finement, or binding of the originally massless gluons into 
massive glueballs. Zeros of the ghost propagator, on the 
contrary effectively supply massless degrees of freedom 
into the equation of state. Moreover, the very existence 
of the mass gap in confining theories requires vanishing 
of the overall contribution of massless degrees of free- 
dom from the low-temperature equation of state. Hence, 
there arises the constraint (|22p on the infrared exponents. 
Confinement resolves itself into a new kind of analyticity. 
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